Microscopic chaos from Brownian motion? 



In a recent Letter in Nature, Gaspard et al. |T|] claimed to present empirical 
evidence for microscopic chaos on a molecular scale from an ingenious exper- 
iment using a time series of the positions of a Brownian particle in a liquid. 
Q\ • The Letter was preceded by a lead article emphasising the fundamental 

0^ . nature of the experiment. In this note we demonstrate that virtually iden- 

_j , tical results can be obtained by analysing a corresponding numerical time 

p series of a particle in a manifestly microscopically nonchaotic system. 

^ I As in Ref. we analyse the position of a single particle colliding with 

many others. We use the Ehrenfest wind-tree model |0] where the point- 
like ( "wind" ) particle moves in a plane colliding with randomly placed fixed 
CN ■ square scatterers ("trees", Fig. la). We choose this model because coUi- 

^ ! sions with the flat sides of the squares do not lead to exponential separation 

t:^- I of corresponding points on initially nearby trajectories. This means there 

O ■ are no positive Lyapunov exponents which are characteristic of microscopic 

Q ■ chaos here. In contrast the Lorentz model used in as an example simi- 

. lar to Brownian motion is a wind-tree model where the squares are replaced 

by hard (circular) disks (cf.[|l|, Fig. 1) and exhibits exponential separation 
of nearby trajectories, leading to a positive Lyapunov exponent and hence 



■ microscopic chaos. 

Q ■ Nevertheless, we now demonstrate that the nonchaotic Ehrenfest model 

^ . reproduces all the results of Ref. A particle trajectory segment shown in 

Q ! Fig. lb is strikingly similar to that for the Brownian particle, (cf.[jl|. Fig. 2). 

Our subsequent analysis parallels that of Ref. [0, where more details may 
be found. Thus the microscopic chaoticity is determined by estimating the 
j_j ■ Kolmogorov- Sinai entropy hj^s, using the method of Procaccia and others 0, 

^ ^ via the information entropy K{n, e, r) obtained from the frequency with 

which the partical retraces part of its (previous) trajectory within a distance 
e for n measurements spaced at a time interval r. Since for the systems 
considered here h^s equals the sum of the positive Lyapunov exponents, 
the determination of a positive hxs would imply microscopic chaos. As 
in we find that K grows linearly with time (Fig Ic and |l[] Fig. 3), giving 
a positive (non-zero) bound on hxs (Fig Id and |jl| Fig. 4). Indeed our 
Figs. Ib-d for a microscopically nonchaotic model are virtually identical with 
the corresponding figures 2-4 of |1]]. Therefore Gaspard et al. did not prove 
microscopic chaos for Brownian motion. 

The algorithm of [HI H] as applied here cannot determine the microscopic 



1 



chaoticity of Brownian motion since the time interval between measurements, 
1/60 s in |l|] , is so much larger than the microscopic time scale determined by 
the inverse collision frequency in a liquid, approximately 10~^^ s. A decisive 
determination of microscopic chaos would involve, it seems at the very least, 
a time interval r of the same order as characteristic microscopic time scales. 
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Figure 1: Brownian motion results of ^ numerically reproduced from the 
nonchaotic Ehrenfest wind-tree model (notation as in The square scat- 
terers have a diagonal of 2 length units and fill half the area considered. The 
particle moves with unit velocity in four possible directions. The position on 
its trajectory is determined for 10^ points separated by one time unit, (a) 
Two nearby trajectories split only at a corner C; no exponential separation 
occurs (cf [H, Fig. 1). (b) A typical trajectory is diffusive, with an cj"^ power 
spectrum (inset), cf [|I|, Fig. 2. (c) The information entropy K{n,e,T) for 
r = 1 and e = 0.316 x 1.21" with m = . . . 25, cf [I, Fig. 3. (d) The 
envelope of the slopes of these i^-curves, h{e, r) appears to imply a positive, 
ie. chaotic, hxs for the Ehrenfest model, as for Brownian motion (cf. |jl|. 
Fig. 4). 



